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Abstract 

The Lemaitre and Schwarzschild analytical solutions for a relativistic 
spherical body of constant density are linked together through the use of 
the Weyl quadratic invariant. The critical radius for gravitational collapse 
of an incompressible fluid is shown to vary continuously from 9/8 of the 
Schwarzschild radius to the Schwarzschild radius itself while the internal 
pressures become locally anisotropic. 
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1 Introduction 



As early as in 1933, G. Lemaitre (IJ emphasized that the most general stress- 
energy tensor associated with a spherical distribution of matter is locally aniso- 
tropic, i.e. expressed in terms of different radial (p r ) and tangential (p t ) pres- 
sures. This generalization which goes beyond the standard perfect fluid ap- 
proximation (p r = Pt) has been extensively used in the recent literature about 
equilibrium Q and collapse || of very compact objects. In particular, the 
Lemaitre case of a collapse with vanishing radial pressure (p r — ; pt ^ 0) 
nicely illustrates how physical (naked) singularities might naturally arise in the 
context of classical general relativity. 

In this note, we will restrict ourselves to the study of equilibrium conditions 
for an incompressible, spherical object sustained by different radial and tangen- 
tial pressures (with p t > p r ). This is what we call the Lemaitre-Schwarzschild 
problem. First, we briefly remind the reader of two analytic solutions (pt = p r 
and p r — 0, respectively) to outline the generic behaviour expected from a dom- 
inant tangential pressure (pt > p r )- Then, we introduce the useful quadratic 
Weyl invariant to link these two extreme solutions and to determine the real 
nature of the so-called Schwarzschild singularity. This relevance of the Weyl 
tensor in the study of gravitational collapse is therefore complementary to pre- 
vious works on conformally flat anisotropic spheres . 

Finally, we illustrate numerically how the boundary of a compact sphere can 
vary continuously from 9/8 of the Schwarzschild radius to the Schwarzschild 
radius itself, without breaking the equilibrium conditions. We confirm there- 
fore that the gravitational redshift of the radiation emitted at the surface of a 
supradense inhomogeneous star could, in principle, be unbounded. 



2 Equilibrium of locally anisotropic spheres 

The interior metric for a static spherically symmetric distribution of matter can 
be written asF] 

ds 2 = e 2 ^dt 2 - e 2X ^dr 2 - rW (1) 

where r is the radial Schwarzschild coordinate and dfl, the solid angle element 
(dn 2 = d0 2 + sin 2 6d<p 2 ). 

As first shown by Lemaitre , the most general spherical distribution of matter 
bounded by gravitation is locally anisotropic. This remarkable feature can be 
directly deduced from the null divergence of the stress-energy tensor 

\p = (2) 
which describes an "anisotropically sustained" bodyQ 

1 We will conventionally take the god-given units where c = 1- 

2 We prefer this expression to the frequently used "pressure anisotropy" which could bring 
the confusion that the pressure would depend on angular variables while, in fact, the radial 
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Tf = diag(p{r), -p r (r), -p t (r), ~Pt{r)), (3) 

with 

T 

Pt-Pr = ^ [Pr + (P + Pr)v'\ , (4) 

the so-called Tolman-Oppenheimer-Volkoff equation. In fact, this relation ([!]) 
directly derives from the radial component of Eq. (Q) but can also be obtained 
from the well-known Einstein equations^] relating the metric ([!]) to the stress- 
energy tensor (|^) (see || for an explicit writing of these three equations). 
Lemaitre only treated the special case of vanishing radial pressure p r in his 
favourite (elliptic) geometry. The full integration of the three independent Ein- 
stein differential equations actually requires state equations for the density and 
pressures, as well as boundary conditions. For recall, the famous exterior metric 

ds 2 = (l - y) dt 2 - (l - y) ~* dr 2 - rW (5) 

with 

rs = 2M, (6) 

the Schwarzschild horizon associated with the sphere of mass M and radius r\ , 
results from the p(r > r%) = p{r > r\) = state equations and g^ir — > oo) = 
r)iiv boundary conditions (G = c = 1). 

Here, the hypothesis of uniform energy density p(r < r\) — po together with a 
regular metric at the origin of coordinates allow us to integrate one among the 
the three differential equations to obtain 



e 



-2A 



K 2 

with 



= 1-55 00 



n' 



n 2 = :i. (8) 

rs 

However, we are still left with two differential Einstein equations involving the 
function v(r), namely Eq.(Q) and 



2 V Po 

Consequently, one additional equation of state for the radial and (or) tangential 
pressures is needed to completely determine the interior metric defined by Eq. ([!]) 
and to analyse the equilibrium conditions of positive and finite pressures 

< p r ,t < oo- (10) 



pressure is just locally different from its tangential counterpart. 

3 Two of the authors (D.L. and A.F.) have shown how to include the effect of the cosmo- 
logical constant as a relique of an underlying bosonic string theory (cf. |5|). But, we will not 
consider that general case here as we have proved that it was not consistent with data of 
precessing pulsars binary systems. 
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2.1 The Schwarzschild perfect fluid (1916) : p r = p t 

If we now assume locally isotropic pressures, Eqs.([|) and (||) can easily be 



integrated to give 



1/2 , N 1/2 

rs 



and 

Mr) = P0 ^ U2 \ 1/2 V 7 1/2 - (12) 

The integration constants are indeed fixed by imposing to recover the exterior 
Schwarzschild solution (see Eq.(^)) at the boundary r = rv The equilibrium 
condition (see Eq.(|l(i|)) for the spherical distribution of a perfect fluid is ob- 
viously violated at r = 0, once the radius reaches a critical value located 
slightly above the Schwarzschild radius r s : 

r? IN = Irs- (13) 

Below this minimal value, the central region of the sphere begins to collapse such 
that the redshift of the radiation emitted at the surface of the sphere cannot 
exceed the critical value of 2. 



2.2 The Lemaitre vaults (1933) : p r = ; p t ^ 

In a similar way, if the material sphere is only sustained by its tangential pres- 
sures, Eqs.(||) and (j^) imply 

e->=(l-^) ! (l-^', (14, 

and 

So, as already stated by Lemaitre in the early thirties, one can build a sphere 
only supported by transverse tensions, at the manner of vaults in architecture. 
At this point, it is historically fair to recall his analytical result |l| 

p t (r) = ^tan 2 ( X ) (16) 

where \ 1S a coordinate used to locate a point on the S 3 space. The trivial 
change of variables r — > 72.sin(%) allows us to convert Eq.(|l5|) into Eq.(|l6|). 
These equations tell us that, here, the tangential pressure never vanishes except 
at the origin and is singular at the boundary when ri tends to rs- Consequently, 
the critical radius for equilibrium is now the Schwarzschild radius itself: 

r^ IN =r s . (17) 
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2.3 Prom Schwarzschild to Lemaitre : p t > p r 

If the tangential pressure is everywhere smaller than the radial one, then the 
boundary condition p r (Xx) = implies Pt(r±) < and thus nothing could pre- 
vent the sphere of radius r± from collapsing. 

On the other hand, if pt > p r , we expect a smooth transition from |rs to r$ 
for the minimal radius, when the radial pressure tends to zero. Our expecta- 
tion is indeed supported by the corresponding smooth evolution of the maximal 
rcdshift factor 

zmax = e-^" JN > - 1 (18) 

of the radiation emitted at the critical surface. As already stated, this critical 
factor is equal to 2 in the locally isotropic fluid approximation (Schwarzschild). 
But Eqs.(p^) and (|l7]) tell us now that this measurable quantity diverges when 
the pressure p r goes to zero (Lemaitre). Would this mean that we only face a 
naked singularity in the extreme case of Lemaitre's vaults? To answer this ques- 
tion, one usually considers the behaviour of the quadratic invariant associated 
with the Riemann tensor, the so-called Kretschmann scalar. 
For the Schwarzschild inner solution, the Kretschmann scalar is given by 



On the other hand, for the Lemaitre solution, we obtain 



(iw^ p<T ) L = 4 ( 5 + 4 - + 12 4 y (*>) 



So, in these two extreme models, a pole in the Kretschmann scalar simply reveals 
the seat of infinite pressure when the boundary of the material sphere reaches the 
minimal radius allowed for equilibrium. It is therefore tempting to conjecture 
that this remarkable feature holds true in all the p t > p r intermediate models. 
The use of the quadratic invariant associated with the Weyl component of the 
Riemann tensor turns out to be a rather powerful tool to prove that a quasi- 
static collapse is indeed triggered when the radius of the sphere lies somewhere 
between frg and r$- 



3 The Weyl quadratic invariant 

With our choice of the metric signature (+, — , — , — ), the Weyl tensor, defined 
as the non-Ricci irreducible component of the Riemann tensor, is given by the 
following formula M: 



1 

d~ 2 

R 



- g^R 



(d-l)(d-2) 



(s w 9v<j - g^gvp) ■ (21) 
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This tensor is traceless and non-trivial if the dimension d of space-time is at 
least equal to four. It is responsible for tidal effects in general relativity but 
vanishes, for example, in all the Friedmann-Lemaitre-Robertson- Walker homo- 
geneous Universes. 

The Weyl quadratic invariant associated with the static metric given in Eq.(|l|) 
reads 

W = W^ pa W^ = I j/cGv - Pt -p)+3 } ( 22 ) 

with k — 8nG- If we deal again with an incompressible fluid of uniform density 
po, Eq.(Q) implies the remarkably simple relation 

W -§(^ (23, 

expressing the fact that the pressure difference p r — pt induces then local tidal 
effects in the anisotropically sustained sphere of matter. 

As it should be, the Weyl quadratic invariant identically vanishes in the 
Schwarzschild homogeneous model [p = constant, pt = p r ) but varies indeed like 
the square of the tangential pressure in the case of Lemaitre's vaults (p r = 0). 
Following our conjecture, we expect therefore a singularity at the place where 
the tangential pressure {pt > p r ) diverges. 

For illustration, let us consider the simple ansatz of Bowers and Liang Q]: 



and 



(i_4y_(i_mV 

consistent with Eqs.(|7|) and (||). The third differential equation given in Eq.(|j) 
implies then 

11 r 2 / r 2 V 1 

Pt-Pr = ^— (g -?)(3p r + Po)Or + Po)^2 ( 1 ~ ^2 ) ' ( 26 ) 



such that the Schwarzschild (Eqs. (|ll| , |l2| )) and Lemaitre (Eqs.( [l4|l5[ )) models 
are easily recovered in the limits q — > ^ and q — > 0, respectively. For the 
intermediate case p t > p r under scruting, a more careful study of Eqs.(|25|) and 
( p6| ) reveals (see, for instance, Fig.l) that instability always occurs first at the 
center of the sphere if < q < 1/2- Indeed, for q ^ 0, the tangential pressure 
diverges at r = when the radius r\ tends to the following critical value 

t min = r s _ (2?) 
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This interpolation of Eqs.(13) and (17) is nicely confirmed by a numerical anal- 
ysis of the Weyl quadratic invariant based on Eqs.(23) and (26) and presented 
in Figj|. Notice that our use of this invariant allows us to conclude that the 
glimpse of a naked singularity can only be caught in the Lemaitre case (q = 0). 
For the other cases (0 < q < 1/2), the collapse always begins at the center 
where the tangential and the radial pressures start to diverge simultaneously. 
In fact, tidal forces cancel out at the center r = (see Fig.||) as expected from 
the Gauss theorem for a spherically symmetric gravitational field. 



4 Gravitational collapse of an incompressible fluid 

In this section we present a simplistic scenario of gravitational collapse to illus- 
trate the important role played by the Weyl quadratic invariant in the physics 
of locally anisotropic spheres. 

The rather smooth behaviour of the Weyl quadratic invariant (see Fig.|] 
and Eqs.(23 2^)) suggests that the gravitational collapse of a Schwarzschild 
supradense body (with q = 1/2) towards the Lemaitre fleeting singularity 
(q = 0) might in principle occur under special circumstances. 
Given the difficulty in deriving realistic equations of state for anisotropic matter, 
the common procedure is to specify an ad hoc relation between the radial and 
tangential pressures. Let us therefore consider the collapse of an incompressible 
fluid and assume that the radial and tangential pressures start to adjust among 
themselves in order to remain both finite. Within such a prescription, we pre- 
clude an early implosion once r\ reaches 9/8 rs (see Fig. 3 for illustration). In 
other words, we impose on the Tolman-Oppenheimer-Volkoff equation (||) to be 
fulfilled at any time and, thus, at any new position taken by the boundary of 
the body. 

At that point, we would like to emphasize that all those assumptions arc rather 
speculative in the sense that we do not propose any particular physical mecha- 
nism that would lead to them (in particular, the density of the fluid could hardly 
remain constant in a real collapse). However, this set of simplistic hypotheses 
will help us in showing the utility of the Weyl quadratic invariant. 
In order to maintain the standard equilibrium condition (0 < p r _t < oo) despite 
the gravitational collapse, the radial pressure p r has to die away while the tan- 
gential pressure decreases near the center and increases near the surface, making 
the sphere looking more and more like a droplet strongly curved by its surface 
tension (as the anisotropy parameter q is lowered to zero). This change of in- 
ternal constitution by local anisotropisation of pressures prevents a premature 
implosion of the sphere. 

The tidal forces appear first far from the surface and grow up as the quasi-static 
collapse goes on. The maximum intensity of those tidal forces moves then from 
the heart of the sphere, when it is still quite isotropically sustained (q ~ 1/2 ; 
r\ w 9/8 rg ; p r w p t ), to its boundary, when the sphere becomes a Lemaitre 
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vault (q « ; r± w rs ; p r — * 0). In this way, the shrinking body would be 
more and more tangentially sustained by (local) tidal forces until the boundary 
eventually reaches its Schwarzschild radius where implosion into a black hole 
necessarily happens (see Fig.Eh. 

A possible astrophysical support in favour of such a scenario involving anisotrop- 
icaly sustained spheres would of course be the direct observation of redshift 
factors 

z max = ( 3 )£ _ j (28) 
much larger than 2, although not related to the expansion of the Universe. 

5 Conclusion 

By analysing the behaviour of the rather powerful Weyl quadratic invariant (see 
Eq.(p2|)), we have illustrated how the gravitational collapse of a spherical body 
down to its Schwarzschild radius is in principle possible through the appear- 
ance of inner tidal forces whose intensities are proportional to the local pressure 
anisotropy. 

Consequently, the Lemaitre's original hypothesis concerning tangentially sus- 
tained bodies could be of interest for further investigations in astrophysics and, 
in particular, for the study of internal structure of supradense bodies from the 
direct observation of sizeable (non-cosmological) redshifts. Needless to say that 
precise state equations have first to be derived in order to proceed. Indeed, 
such bodies would rather be constituted of degenerate nuclear matter with lo- 
cally anisotropic pressures (see ||). 

The careful reader has certainly noticed that the analytical expression given 
in Eq.(|2^) for the Weyl quadratic invariant is also valid in the dynamical case 
where both the metric and matter fields depend on time t and radial coordinate 
r. Further work on gravitational collapse of locally anisotropic spherical bodies 
in terms of poles in this invariant would therefore be worth pursuing. In partic- 
ular, the regularity of the Weyl quadratic invariant can be used to characterize 
the nonsingular asymptotically flat solutions to the static spherically symmetric 
Einstein- Yang-Mills equations with SU{2) gauge group. 
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Fi glire 1 '. Radial and Tangential Pressures for different values of the anisotropy parameter 
q (ri = 1.3r s ; r s = 1). 
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Figure 2: Weyl Quadratic Invariant W (see Eq. (|2^ ), related to the intensity of tidal effects, 
for different values of the anisotropy parameter q. It is proportional to the square difference 
of the two preceding figures (r± = 1.3rg ; rg = 1). 
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Figure 3: Modification of the radial and tangentiai pressures (p r , pt, respectively) needed 
when the limit r\ of the sphere goes below the critical radius |rg of the isotropic case where 
Pt = Pr- (rs = 1 hi the figure.) 
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Figure 4: Weyl Quadratic Invariant (W.Q.I.) given by Eq. ( [2j ) as a function of the anisotropy 
parameter q and the position r when the limit n of the sphere goes below the critical radius 
|rs of the isotropic case (pt = p r ). (rs = 1 in the figure. J 
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